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ABSTRACT 



Aims. We investigate the mass loss of highly evolved, low- and intermediate mass stars and stellar samples with subsolar metallicity. 

We give a qualitative as well as quantitative description which can be applied to LMC/SMC-type stellar populations. 

Methods. For that purpose we apply the same approach as we did for solar metallicity stars and calculate hydrodynamical wind 

models including dust formation with LMC and SMC abundances under consideration of an adapted model assumption. In particular, 

we improved the treatment of the radiative transfer problem in order to accommodate larger non-local contributions occurring with 

smaller opacities. For each wind model we determine an averaged mass-loss rate. The resulting, approximate mass-loss formulae are 

then applied to well-tested and calibrated stellar evolution calculations in order to quantify the stellar mass loss. 

Results. The dynamical models for LMC and SMC metallicity result in mass-loss rates of the same order of magnitude as the solar 

metallicity models which is in this basic approach in agreement with observations. The hydrodynamical properties like e.g. the outflow 

velocity differ (for fixed C/0 abundance ratio) noticeably, though. While critical luminosities of LMC and solar metallicity models 

fairly coincide, the SMC models need higher luminosities to develop dust-driven winds. 

Key words. Hydrodynamics ~ Stars: AGB and post-AGB - Stars: carbon - Stars: evolution - Stars: late-type - Stars: mass-loss 



1. Introduction 

On the asymptotic giant branch (AGB) stars with low- and in- 
termediate main sequence mass (~1 ... 8 Mq, precise values 
depend on the initial element abundances / metallicity and are as 
well model dependent) are in their late stage of evolution. This 
phase is characterised by instabilities leading to stellar pulsa- 
tions as observed in long-period variables (LPVs) or Mira stars. 
ICioni et aL I (120011) . for instance, found that 65% of the stars in 
the considered sample of giants in the Large Magellanic Cloud 
(LMC) are LPVs. Furthermore, this phase is connected to strong 
mass-loss of more than a few 10"^ Mq yr"'. By means of their 
mass loss, these stars contribute significantly to the enrichment 
of the interstellar matter with material processed through nu- 
clear burning reactions in the form of molecules and dust, and 
are therefore not only interesting as objects of stellar evolution 
but also in terms of galactic chemical evolution. In this paper we 
extend our set of models with solar element composition (apart 
from the C/O ratio) to subsolar metallicities. We are using abun- 
dances as observed in the Magellanic Clouds (see table A.l in 
[Helling et al., 2000). The metallicity of the Large Magellanic 
Cloud (LMC) is Z = 0.008, i.e. roughly a good third of the solar 
value, while the Small Magellanic Cloud (SMC) has on average 
an even lower metallicity of Z = 0.004. 
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2. Dust-driven Wind IVIodels 

We base our description of the mass loss of late AGB stars with 
less than solar metallicities on time dependent hydrodynamical 
wind models that include the formation of dust grains which, 
by radiation pressure, drive the massive o utflows. These model s 
have been developed in the last decade by Gauger et al. ( 199C|, 
FleischeretaL (.1992) . .Winters et al., (,1994., ,200Q) . .Jeong et al] 
(I2003h . 

For (apart from the C/O abundance ratio) solar abundances, 
these models yield time averaged mass loss rates of up to a few 
10"^ Mo yr"', both for oxygen-rich and carbon-rich mixtures. 
This is consistent with rates inferred from, e.g. CO rotational 
line observations of AGB stars in the solar neighbourhood. 

The models are obtained from the solution of the non- 
linearly coupled system of equations describing the hydrody- 
namical and thermodynamic al structure of a spherically sym- 
metric, pulsating circumstellar shell, its chemical composition, 
and the nucleation, growth, and evaporation of dust grains. The 
hydrodynamical structure (mass density p and outflow velocity 
v) follows from the equation of continuity and the equation of 
motion which includes radiation pressure on dust grains. The 
law of energy conservation and radiative transfer determine the 
temperature structure. In this work we consider models with 
carbon-rich composition assuming that oxygen is locked in the 
CO molecule and assuming chemical equilibrium. The forma- 
tion, growth, and evaporation of carbon grains is calculated ac- 



2 



A. Wachter et al.: Dust-driven winds and mass loss of C-rich AGB stars with subsolar metallicities 



Table 1. Essential input parameter ranges covered by our collec- 
tion of subsolar metallicity models. 



type 


LMC 


SMC 


[Le] 


5000. 


. . 15000 


7000. 


. . 12000 


T. [K] 


2200. 


. . 3200 


2400. 


. . 3000 


M, [Me] 


0.7 .. . 


1.0 


0.6. . 


1.1 


fc/eo 


1.3/1. 


8 


1.8 




P[d] 


325.. 


.650 


450.. 


.650 


Av [km/s] 


2/5 




5 





cordin g to the moment method developed by iGail & Sedlmayd 
(119881) and'GaugeretaD 11990^. 

The equations are discretised in Lagrangian coordinates, i.e. 
in a co-moving coordinate system. A detailed description of 
the tra nsfor mation an d discre tisation is given in [Fleischer et alJ 
(119921) and iFleisched (11994 . Furthermore, a summary of the 
latest p hysical assumptions can be found in ISchirrmacher et al] 



When applied to winds with subsol ar metallicities thes e 
models show a rather quasi-static nature ( iHelUng et all |2002[) . 
As already pointed out by these authors, this is not consistent 
with observations of AGB stars, e.g. in the Magellanic Clouds 
which indicate mass-loss rates of the same order of magnitude 
as they are observed from stars with solar element abundances. 
Therefore, we reconsidered the model assumptions. In particu- 
lar, we improved the treatment of the radiative transfer as fol- 
lows: 

In the solar metallicity case the shell has been considered 
to be optically thick so that the flux-averaged dust extinction 
is determined by the Planck-mean absorption coefficient at the 
local equilibrium temperature. However, if the atmosphere is 
considered to be partly transparent, the dust grains are to some 
extent exposed to the direct radiation from the stellar photo- 
sphere. Therefore, the flux-mean dust opacity for lower metallic- 
ity is represented by assuming a «o«-local radiation field, which 
is characterised by a temperature T determined according to 
j4 _ T^e^'^ + T^^{1 - e"'"), i.e. the radiation field is interpo- 
lated in optical depth t between the local equilibrium radiation 
and the stellar (black-body) radiation field (represented by Tgq 
and Ti,, respectively). 

The gas opacity Kg/p is set to a constant value of 2 x 
10~"^ c m^ g ' in the solar models. From figure 2 of lHelling et alj 
( I2OOOI) it is clear that the mean gas opacity drops with lower 
metallicity. Therefore, we reduced the value of the gas opacity 
in the LMC and SMC models to 1 x 1 0"**, and 0.5 x 1 0"'* cm^ g" ' , 
respectively. 

3. Results 

Since the main purpose of this work is to derive mass-loss rates 
for subsolar metallicities, the focus when choosing the model 
parameters was primarily on covering the range of stellar pa- 
rameters M, L, and T provided by the stellar evolution code for 
each fixed element composition (LMC/SMC) than on perform- 
ing a detailed parameter study. Nevertheless, the dependence on 
the carbon overabundance is shortly discussed in the following. 
Table [1] summarises the parameter ranges of our subsolar metal- 
licity model sets. 

Concerning the radial structure of our wind models there ba- 
sically exist two different types of models: 

A) models where effective dust formation and growth takes 
place. In this class of models the radiative acceleration on dust 



exceeds the local gravitation in the region close to the stellar 
photosphere (typically inside of » 5/?*). In this case, the matter 
is lifted out of the stars' gravitational field by radiation pressure 
on dust, and a layered structure of the hydrodynamical and the 
dust quantities results in the inner part of the circumstellar shell 
(see, e.g., figure [T]|. Generally this type of models shows high 
mass-loss rates (> a few lO^^Moyr"') and outflow velocities in 
excess of » 5 km s" ' . 

B) models showing a smooth, quasi-stationary wind struc- 
ture. In these model, radiation pressure on dust is not sufficient 
to lift the matter out of the stars' gravitational field, the resulting 
mass loss rates (and outflow velocities) are significantly smaller 
than in the type A models. 

Th is separation has been discussed in detail in [Winters et al.l 
(I2000h for the case of solar element abundances. 

As the distinction between A-type and B-type models is due 
to the ratio of radiative to gravitational acceleration (larger or 
smaller than unity) in the inner shell region, the same type of 
separation is found for sub-solar metallicities. For the SMC and 
LMC models, this separation is later on considered in terms of a 
"critical luminosity". 

In order to characterise the outflow of each time dependent 
wind model time-averages of several quantities are calculated. 
The most obvious quantities for this purpose are the mass-loss 
rate (M) and the final outflow velocity (voo). Furthermore we 
consider the dust-to-gas ratio (pd/Pg) and the ratio of radiative 
to gravitational acceleration {a}, since these are expected to be 
significant for the general behaviour of the wind. 

To calculate the averages the same procedure as used in the 
solar metallicity models is followed, i.e. we first average M, v, 
and pd/g in the radial coordinate between 40 and 60 and then 
over time, typically 20 periods for LMC and 80 periods for SMC 
models. For {a} a different procedure is followed, though. It is 
a time-average only, namely the average of the ratio of radiative 
to gravitational acceleration at the radius of the first maximum 
of the condensation degree, i.e. in the innermost dust layer (see, 
e.g., the innermost peak of a in the second panel of figurelTJ. 



3.1. Influence of the CIO ratio 

The C/O abundance ratio determines how much carbon is avail- 
able for dust formation, since we consider carbon rich chem- 
istry and total CO-blocking. As far as the derivation of a mass 
loss formula applicable for stellar evolution is concerned, the 
solar metallicity models already demonstrated that the mass 
loss rate does not depend heavily on th is model parameter (see 
lArndt et al.L IT997I: IWachter et all l2002h . Instead it can be seen 
as a critical parameter in the sense that the mass loss rate re- 
mains the same order of magnitude once the C/O ratio is above 
a certain threshold. 

The results of a set of calculations for LMC and SMC abun- 
dances with varying C/O ratio are summarised in table |2] In the 
LMC models with ec/eo = 1.3 and 1.5 the material expands so 
slowly that after 110 and 90 periods the outermost grid point 
has just reached 10 and 25 stellar radii, respectively. This is well 
below the region of 40 to 60 stellar radii where we generally 
take the averages. The calculations have not been followed any 
further since the models up to that instant show a virtually sta- 
tionary radial structure. 
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Table 2. Variation of the carbon overabundance ec/o for other- 
wise fixed parameters. 

fc/o I periods \M\ (a) v„ Pd/g 



type: LMC 


M = 0.8 Mo, r = 2800 K, L = 15000 Lq, 
P = 400 d, Av = 5 km s"' 


1.30 
1.50 


Expansion after 1 10 P: < 10 R* 
Expansion after 90 P: < 25 Ri, 


1.70 
1.80 
2.50 


70-90 3.51e-05 6.10 18.83 2.01e-3 
70-90 3.30e-05 3.92 21.04 2.68e-3 
70-90 4.58e-05 6.03 30.16 7.13e-3 


type: SMC 


M = l.OMe,T = 2600 K, L = 1 1000 Lq, 
P = 650 d, Av = 5 km s"' 


1.30 

a < 1 1.50 
1.70 


Expansion after 150 P: < 35 Ki, 
70-150 7.17e-07 0.06 3.23 1.69e-4 
70-150 4.74e-06 0.11 4.95 1.17e-3 


2.00 


70-150 3.69e-05 1.26 8.56 1.82e-3 
70-150 4.55e-05 2.19 10.61 2.75e-3 


type: SMC 


M = 0.8 Mg, r = 2800 K, L = 10000 Lq, 
P = 600 d, Av = 5 km s"' 


^^r 1-30 
1.50 


130-150 1.24e-07 0.27 1.91 1.93e-9 
70-150 8.90e-07 0.21 3.75 4.74e-4 


1.70 

a>l 1.80 
2.00 


70-150 1.72e-05 1.02 8.31 1.27e-3 
70-150 2.88e-05 2.04 10.57 1.75e-3 
70-150 3.76e-05 3.11 12.10 2.89e-3 



Table 3. Mass-loss rates of SMC, LMC, and solar wind models 
with identical input parameters; ec/o - 1.8, Av = 5 km s"', P 
= 500 d at 8000 La, 600 d at 10000 Lg, and 700 d at 12000 Lq. 
Mass loss averages are over time interval 70-150 periods, or (1) 
70-90, (2) 70-210. 



M 


T, 






|M|[10-^Moyr- 


-'] 


[Mo] 


[K] 


[Le] 


SMC 


LMC 


solar 


0.8 


2400 


10000 


2.19 


1.03 


5.45 


0.8 


2600 


8000 


2.75 


3.31 <" 


2.38 


0.8 


2600 


10000 


5.14 


8.18 <" 


3.57 


0.8 


2800 


10000 


2.88 


4.58 


2.33 


0.8 


2800 


12000 


4.69 


6.74 


3.23 


0.9 


2800 


10000 


1.92 


3.17 <" 


1.94 


1.0 


2600 


10000 


1.19 


3.81 <" 


2.190 



3.2. Comparison of models with different element 
abundances 

Generally, when comparing the Magellanic Cloud (MC) models 
to the solar ones with identical input parameters one finds lower 
final outflow velocities and lower dust-to-gas ratios. The radia- 
tive (in terms of the gravitational) acceleration is weaker as well. 
For the mass-loss rate, there are input parameters where the MC 
models show higher values than the corresponding solar ones, as 
can be seen in table[3] The fact that one finds higher outflow ve- 
locities with increasing heavy element abundance, and mass-loss 
rates of the same order of magnitude is not surprising consider- 
ing results from basic stellar wind theory. The same behaviour is 
seen assuming an isothermal wind with an outward directed r"^ 
force applied outside the critical point. There the mass-loss rate 
is set as long as the force is only acting in the supersonic regime, 
that is outside the critical point. 

Typically the outflow velocities of the solar metallicity mod- 
els are higher by about a factor of 2.2(+0.2) than those of the 
LMC, and 4(+l) than those of the SMC. The dust-to-gas ra- 
tios are higher by about a factor of 1.3(±0.1) and 2.3(+0.2), re- 



spectively. For the acceleration the trend is less clear, i.e. the 
data show a wider spread, especially in the SMC case. The 
solar abundance models produce values higher by a factor of 
3.6(+ 1 .2), and 5.8(+ 1 .2), respectively. 

Figure [T] depicts the radial structure of the model with M = 
0.9 Mq, = 2800 K, L = 10"* Lq, P = 600 d, Av = 5 km s"', 
and ec/eo - 1.8 for SMC, LMC, and solar abundances at the 
instant of 90 periods after starting the model calculation. It gives 
a good representation of these trends, especially for the outflow 
velocity and the dust density, recalling that they are not only 
time but as well radial averages. The acceleration trend is more 
indirectly represented, since the graph is a time snapshot, but the 
acceleration average is a time average of the value at the radius 
of the first maximum of the condensation degree. 

In the lowest panels two more trends can be seen. The con- 
densation degree which is closely related to the dust-to-gas ra- 
tio by mass according to Pd/Pg ^o(^ - l)/c (see eq. (4.7) 
of IWintersl 1 1 9941) is less pronounced and more irregular in LMC 
and SMC models than in solar models. While in this solar metal- 
licity model the condensation degree is virtually 1 in the dust 
shells, the SMC model only reaches maximum value of about 
0.8. That is all the condensible material is accumulated in dust 
grains in the solar model, but with lower metallicity the con- 
densation is less efficient. The second trend resulting from this 
simple graphical comparison of respective models concerns the 
average grain size. In the SMC case the grains reach much larger 
radii than with solar abundances. Since the nucleation rate drops 
with decreasing metal abundance, there are fewer grains in that 
case to collect the available C atoms. Additionally, the grains 
have more time for growing due to their lower velocity, hence 
they become larger. 

3.3. Approximative mass-loss formulae 

To give a mass-loss description applicable to stellar evolution we 
restrict our model set to those of type A, describing stable dust- 
driven wind. These are characterised by the criterion (a), > 1 . 

Furthermore, one has to consider the treatment of the "me- 
chanical" input parameters Av and P as well as the C/O abun- 
dance ratio for this purpose. Fortunately, the models are most 
sensitive to the stellar parameters luminosity L, effective tem- 
perature T, and mass M (for the so lar metallicity models already 
pointed out by lArndt et all Il997h . Therefore, we followed for 
the LMC and SMC models th e same approach as for solar metal- 
licity models (I Wachter et al. l 2002) and fixed the velocity ampli- 
tude Av to a value of 5 km s"' . Setting the abundance ratio ec/o 
to 1 .8 allows for the fact that it has to be high enough to produce 
a stable wind without being unreasonable. 

Concerning the period P we relied on observed period- 
luminosity relations of LPVs. Wachter et al. (2002) took the de- 
pendence of the mass-loss rate on the period into account by 
transferring the slight P dependence to an additional luminos- 
ity term. Now, we take the period dependence more directly into 
account by choosing only those models where the P-L relation 
(Groenewegen & Whitelock, 1996) is fulfilled. This is possible 
since we now have a larger set of models at hand. 

For the selected LMC models a multi-linear least square fit 
of the form log M = gq + fi * log xi + . . . was performed leading 
to the following formula: 

|M[Meyr-']| = 3.80 X 10"-"^ X (M*[Mo])"^-^* 
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Fig. 1. SMC (left), LMC (middle), and solar (right) model with identical input parameters M = 0.9 Mq, = 2800 K, L = IO'^Lq, 
P = 600 d, Av = 5 km s ec/eo = 1-8. Shown is the radial structure at time t - 90 P, being the stellar radius at the start 
of the calculation t - QP. The depicted quantities are (left scales) the radial velocity v (solid blue line), the speed of sound cj 
(dotted blue lines - in-/outwards, zero line as guide to the eye), the density p (solid magenta), the nucleation rate per hydrogen atom 
J*ln(H) (solid black), and the degree of condensation of the dust (solid red), (right scales) the gas temperature (dashed red) and 
equilibrium temperature (solid cyan) T, the ratio of radiative to gravitational acceleration a of the material (dashed blue), the dust 
density nd (dashed cyan), and the mean dust particle radius (r) (dashed blue). 



with a correlation coefficient of K — 0.98. 
For the SMC models the result is: 



|M[Moyr- 



-3.01 



2.34X 10"^ x(M*[M0]) 



2600 



104 



(2) 



with a correlation coefficient of A' = 0.97. 

The formula for solar element composition as published in 
IWachter et al.l (|2002[) reads, for comparison: 



|M[Moyr"']| = 4.52 x 10"^ x {M^{Mq\) 



-1.95 



2600 



104 



(3) 



where the correlation coefficient - 0.97. 

As already pointed out for the solar metallicity case, the most 
influential parameter in the formulae for LMC/SMC abundances 
is again the effective temperature. This reflects the extreme sen- 
sitivity of the dust nucleation with respect to the local tempera- 
ture. 

Graphical illustrations of the formulae for LMC and SMC 
are given in figure |2] for a fixed mass of = 1.0 Mq, tem- 
perature range of 2400 . . . 3000 K, and luminosity of 5000 . . . 
15000 Lq. The lines on the base of the diagrams are projections 
of contour lines, i.e. they mark those temperature and luminosity 



values for which the resulting mass-loss rate is lO^^^ 5 x 10"^, 
10"^, and 5 x 10"^ M© yr"', respectively. These show that for the 
SMC higher luminosities / lower temperatures are necessary to 
reach the same mass-loss rate as in the LMC. 

To give the validity limit of the above formulae the selection 
criterion needs to be expressed in terms of the stellar parameters, 
since (a) is a quantity which does not appear in stellar evolution 
calculations. By considering the ratio LjM over temperature T 
for each metallicity set it turns out that the criterion can be ap- 
proximated using oc T^. We find that equations ([TJ and Q are 
valid for 



L > La 



' 5.58 X 10-^ r^ i''** M, LMC 
5.47 X 10-4 M, SMC ^ ' 



The need to include such precise values of the temperature expo- 
nent expresses again the steep dependence of the dust nucleation 
rate on the local kinetic gas temperature. 



4. Results from application to stellar evolution 

In order to get some insight into the overall mass loss of subsolar 
metallicity stars on the AGB compared their solar counterparts, 
we applied the derived mass-loss rates to stellar evolution cal- 
culations. For that purpose, we used again the well-calibrated 
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(a) LMC: 




Fig. 2. Graphical illustration of the fit formulae for the mass loss 
rate for fixed mass = 1.0 M©. 



Cambridge evolution code (cf. IWachter et all l2002l and refer- 
ences therein), as we did for the solar metallicity models. A par- 
ticular advantage of these evolution models is a good match of 
their effective temperatures for stars high on the AGB, accordin g 
to a semi-empirical adjustment made by Schroder et al.l (Il999h . 
This is an important point: Only a small mismatch in Tsg can 
result in a significantly different mass-loss rate, despite the use 
of the very same mass-loss formula, because of the high depen- 
dence of the mass-loss on T^ff of the stellar models. 

Mass loss is included by reducing the stellar mass at the outer 
boundary, whereby the particular formulation depends on the 
current evolutionary stage of the model. Mass loss is considered 
only after the star runs out of hydrogen in the core and reaches 
the red giant branch (RGB). The formulae derived in section [331 
are valid above the respective critical luminosities. These are 
reached when the star has developed up to the tip of the AGB 
and has entered the thermally pulsing phase. Before this phase 
we apply a formula given by Schroder & Cuntz (2005) who have 
reconsidered the well-known Reimers formula dReimersl Il975h 
by using theoretically motivated arguments. As a result it now 
includes two new terms containing a dependence on the effec- 
tive temperature and the surface gravity, respectively. Around 
Lcrit a short transition zone is introduced where the code inter- 
polates between the modified Reimers and the formula based on 
dust-driven wind models to avoid artificial jumps in the mass- 
loss rate. 

We have generated model grids with a metallicity of Z = 
0.01 and Z = 0.001 to represent LMC and SMC stars, respec- 



tively, the closest values for which appropriate necessary input 
data is available. 

Table |4] lists the initial masses and the total mass lost over 
the whole life of the star, according to a number of our evolution 
models with mass loss for subsolar metallicity. Within each grid 
we calculated evolutionary tracks without and with overshooting 
(w ith a value of ^py = 0. 12 for the overshoot parameter as tested 
bv lSchroder etailfT997l) . 

4. 1 . Mass loss histories 

The temporal development of the stellar mass loss rate is of par- 
ticular interest, since the density structure of a planetary neb- 
ula is determined by the mass-loss history of its progenitor 
star Figure [3] shows four examples of such mass-loss histories. 
Depicted are the rates of Z = 0.01 models with an initial mass 
of (a) Mi = 1.25 Mq and (c) Mi = 2.50 Mq, and Z = 0.001 Mg 
models with (b) Mi = 1.95 Mq and (d) Mi = 2.50 Mq. Each 
graph is plotted over the final 200 000 years of the stellar life, 
but the scales of the mass-loss rates differ! The graphs (a) and 
(b) show the results for the lowest initial mass where the criti- 
cal luminosity is reached and a brief onset of a dust-driven wind 
results, respectively, while graphs (c) and (d) show the highest 
initial mass where the evolution is followed until the star leaves 
the AGB. 

What catches the eye at first glance is the fact that the num- 
ber of thermal pulses (seen as spikes) increases and the inter- 
pulse duration decreases with decreasing metallicity (cf. fig- 
ures[3lc) and (d)). This finding is in accordance with the predic- 
tions of other stellar evolution calculations using different mass 
loss prescriptions, see for example Marigo (2001), and refer- 
ences therein. As pointed out there, the exact number of thermal 
pulses depends on the chosen mass-loss formalism. 

Another clear difference is the overall shape of the mass- 
loss histories of models with different metallicity: The Z = 0.01 
graph resembles the mass-los s rates of our models with solar 
metallicity, see figures 2-5 in IWachter et al.l (|2002|) . For these 
metallicities the critical luminosities are fairly similar and the 
stars reside in a regime above that value while they are in the 
thermally pulsing AGB phase. As a result the mass loss in this 
late evolutionary stage is controlled by our formula based on 
dust-driven winds. By contrast, the Z - 0.001, M = 2.5 Mq 
rate is mostly dominated by not dust-driven mass-loss, here de- 
scribed by the modified Reimers law. Only during the last few 
thermal pulses does the star reach luminosities high enough for 
the dust-driven wind formula to apply. Our present models for 
SMC-like metallicity therefore result in a rather smooth and 
broad density distribution of the blown-out material develop- 
ing to a planetary nebul^ Models with solar and LMC abun- 
dances though have distinctly higher mass-loss rates during the 
last ~25 000 years leading to a denser and narrower mass distri- 
bution. 

4.2. Total masses lost 

The total mass lost by a star in the course of its evolution, for 
given initial stellar mass, is an important quantity in modelling 
the galactic chemical evolution. A further quantity of interest 
is the mass of the resulting white dwarf (WD) or, rather, the 
WD mass distribution, since WDs are considered to be the end- 
products of stellar evolution for intermediate masses. The final 

' Naturally, this is a simplified picture, since it disregards different 
velocities with which the material might be driven away. 
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Table 4. Stellar evolution grids: final masses and ages 







Z=0.01 






Z=0.001 






Oov 


- Q 


A 

Oov 


— u. iz 


Oov 


- 


Oov 


— u. iz 




M 


age 




age 


M 


age 


Me 


age 


L^oJ 


LMoJ 


[iU yrj 




[iu yrj 


r A/f 1 

L^oJ 


[10 yij 




[iu yrj 
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stellar (WD) mass is basically identical with the core mass of the 
AGB star supposing that the whole envelope is ejected via stellar 
winds. Therefore we present in the following the total (integrated 
over the whole stellar life) masses lost by stars with subsolar 
metallicity, according to our models. 

Figure m depicts the amount of mass lost by a star of given 
initial mass Mi during its life for the two grids of tracks with 
metallicity (a) Z = 0.01 and (b) Z = 0.001. In addition, we de- 
tail the time-integrated mass-loss of the RGB and AGB phases. 
In both graphs the respective masses lost by models with solar 
metallicity (and same initial mass) are plotted as well, for com- 
parison. It can be seen in figureUa) that the total amount of mass 
loss is virtually the same, but not lower for LMC-like than for so- 
lar abundances. However, the contributions lost on the RGB and 
AGB are shifted, especially for the lower initial masses shown 
more mass is already lost in the RGB phase. For Z - 0.001 fig- 
ure |4jb) shows that the total amount of mass lost is almost the 
same as for solar metallicity for the lower end of initial masses. 
With higher initial mass though, less material is blown away for 
SMC-like abundances. Again there is a shift towards higher con- 
tribution of the RGB, but less pronounced than in case (a). 

5. Discussion 

In order to account for carbon-rich chemistry, we enhance the 
carbon element abundance considering the C/O ratio as a model 
parameter. In this work we considered mainly values of 1.3 



and 1.8, the latter especially in models with SMC abundances. 
This value might be r egarded as rather extreme. Nevertheless, 
Matsuur a et al.l (l2005b reported possible C/O values of ~1.4 for 
LMC stars and even higher for the SMC, resulting from a chem- 
ical equilibrium model to explain the strength of observed C2H2 
absorption features. 

IVan Loon et alj (l2005b observed AGB stars in the Large 
Magellanic Cloud, C-stars as well as M-type stars, and derived 
mass-loss rates of log M[Mq yr"^] ~ -5.6 ... -5 for C-type and ~ 
-5 . . . -4.2 for M-type giants. Furthermore, they give an empirical 
mass-loss formula for these stars dependent on effective temper- 
ature and luminosity and compare their formula to equation (O 
derived from our models with solar abundances. Even though 
the dependence of our LMC formula on effective temperature 
is stronger, it is consistent with their observed value within the 
errors. Furthermore, they argued that even if their observed lu- 
minosity dependence is much lower than our theoretical value, 
they may still be consistent as the effect of decreasing mass-loss 
rate with increasing stellar mass counteracts the luminosity de- 
pendence. 

An issue which remains uncovered by the models presented 
here is related to the assumption of gre y radiative transfer and 
its effect on the resulting mass-loss rates. lHofner et al.l (l2003h in- 
cluded frequency-dependent treatment in their models and found 
those to have cooler and denser upper layers which is in favour 
of a more efficient dust condensation. A radiative pressure 2- 
3 times larger than the grey values results, and the mass-loss 
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Fig. 3. Mass loss histories of stars with different initial masses and metallicities. Labels denote the current mass. 



rates and wind velocities are larger in the non-grey case for the 
three models compared to grey on es calculated with Planck- 
mean opacities. On the other hand, iHelUng et al.l ( |2000 ) inves- 
tigated the effects of molecular opacities on the characteristics 
of carbon-rich dust driven wind models in the frame work of 
grey opacities. These authors compared opacities represented by 
the same constant value used here, with Rosseland-mean and 
Planck-mean averages. As a general result it is found that the 
mass loss rates and wind velocities are highest for the constant 
value of the opacity also used here and lowest for the Planck- 
mean opacities, whereas the Rosseland-mean leads to interme- 
diate values o f the mass loss rate. Considering these results and 
the results of iHofner et alj (|2003|) it is to be expected that the 
inclusion of non-grey opacities in our models may lead to some- 
what reduced mass-loss rates and wind velocities. A quantitative 
assessment of the effect would however require the calculation of 
an appropriate set of models with non-grey treatment of the gas 
opacities for the sub-solar element abundances discussed here. 

Concerning the comparison of h ydrodynamical wind m odels 
with different element compositions iMattsson et afl (l2008l) sug- 
gest that one should rather keep the abundance difference C-O 
fixed than the ratio C/O, since for dust formation the amount of 
free carbon is important assuming complete CO blocking. Doing 
so effectively leads to higher C/O ratios for lower metallicity 



where the oxygen abundance is decreased as well. For the mass- 
loss rates this different approach of comparing models conse- 
quently has not much influence, since the models with higher 
C/O are even further beyond the threshold. For example, we cal- 
culated the LMC and SMC model shown in figure [1] with fixed 
C-O (logio(ec - eo) ~8.72fl instead C/O. This leads to averaged 
mass-loss rates reduced by a factor of 0.7 and raised by a factor 
1.3 (to the value with fixed C/O), respectively. Considering the 
hydrodynamical structure itself, the models with fixed C-O are 
very similar to their solar counterparts, e.g. they show velocities 
of the same order and the degree of condensation reaches 100% 
in the dust shells even in the SMC model. 

Despite the differences in the assumptions of gas opacities 
and radiative transfer we co me to the same conclusion about 
mass loss as Matt sson et al.l ("^OOS): Lower metallicities do not 
necessarily result in lower mass-loss rates. 

As stated before the wind models predict higher critical lu- 
minosities for a dust-driven wind to develop for SMC abun- 
dances than for the LMC, e.g. a factor of 1.7 at log T = 3.5. On 
the other hand, stellar evolution tracks of stars with Z = 0.001 
reach higher luminosities (up to 0.2 dex) at the tip of the AGB 



which corresponds to values of C/O ~ 3.4 and 5.9, respectively 
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(a) Models with metallicity Z = 0.01: 
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(b) ...and Z = 0.001: 
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Fig. 4. Integrated mass loss as a function of initial mass. For each 
initial mass the total mass loss is shown (blue diamonds) as well 
as the respective fractions lost on the RGB (cyan) and the AGB 
(grey). For reference the respective solar values are given as well 
in both plots (red, green, and black crosses). 



than their Z = 0.01 counterparts with same initial mass. In other 
words our models predict SMC stars to undergo mass loss due to 
dust-driven winds, at least those with initial masses of about two 
solar masses and higher, despite the higher critical luminosity 
required to do so. 



6. Conclusions 

In this article we present dust-driven wind models with subso- 
lar metallicities as observed in the Magellanic Clouds. These are 
calculated using an adapted version of the hydrodynamical code 
developed at the Technical University Berlin, with a radiative 
transport description improved for lower opacities. Based on a 
grid of models for SMC and LMC metallicity respectively, we 
derived approximative mass-loss formulae valid above a critical 
luminosity, that is on the tip of the AGB. In order to investigate 
the overall stellar mass loss we applied these descriptions to stel- 
lar evolution using the Cambridge code. 

Our main results can be summarised as follows: 



1. A comparison of the subsolar metallicity models to the cor- 
responding solar ones with fixed C/O ratio reveals that typi- 
cally the outflow velocities of the solar models are higher by 
about a factor of 2.2 than those of the LMC, and a factor of 
4 higher than those of the SMC. This is because in the so- 
lar metallicity case a larger amount of dust is formed which 
can more efficiently accelerate the wind outside of the sonic 
region. 

2. The averaged mass-loss rates of these models are of the same 
order of magnitude as the solar ones. This can be understood 
by the threshold character of the radiative acceleration in a 
dust-driven wind: Only a moderate amount of dust is needed 
to form in the sub-sonic wind region to accelerate the wind to 
super-sonic velocities (equivalent to a > 1 at sub-sonic ve- 
locities). The mass-loss rate is in this case determined in the 
sub-sonic region, irrespective of how much additional dust is 
still formed further out in the wind. 

3. The most influential parameter in the mass-loss formulae is 
the temperature hinting at the sensitive dependence of the 
dust formation to this quantity. 

4. Once the abundance ratio C/O is above a lower limit it still 
has an impact on the hydrodynamical structure of the mod- 
els, e.g. different outflow velocities, but the averaged mass- 
loss rate does not change dramatically anymore. Again this 
is plausible due to characteristics of the driving force in a 
radiation driven wind (see also section l3T2l i. 

5. While LMC and solar abundance wind models result in fairly 
equal critical luminosities necessary to drive a wind by radi- 
ation pressure on dust, higher values are required in the SMC 
models. Since stellar evolution models with lower metallic- 
ity reach higher luminosities on the tip of the AGB, some of 
these stars (with initial masses Mi > 2Mq) suffer dust-driven 
mass loss nevertheless. 
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